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ABSTRACT 

The base of the convection zone is a source of acoustic glitches in the asteroseismic 
frequency spectra of solar-like oscillators, allowing one to precisely measure the acoustic 
depth to the feature. We examine the sensitivity of the depth of the convection zone to 
mass, stellar abundances, and input physics, and in particular, the use of a measurement 
of the acoustic depth to the CZ as an atmosphere-independent, absolute measure of 
stellar metallicities. We find that for low mass stars on the main sequence with OAMq < 
M < I.6M0, the acoustic depth to the base of the convection zone, normalized by the 
acoustic depth to the center of the star, Tcz^n, is both a strong function of mass, and 
varies at the 0.5-1% per 0.1 dex level in [Z/X], and is therefore also a sensitive probe 
of the composition. We estimate the theoretical uncertainties in the stellar models, 
and show that combined with reasonable observational uncertainties, we can expect 
measure the the metallicity to within 0.15 - 0.3 dex for solar-like stars. We discuss 
the applications of this work to rotational mixing, particularly in the context of the 
observed mid F star Li dip, and to distguishing between different mixtures of heavy 
elements. 

Subject headings: stars: abundances, stars: interiors, stars: oscillations 



1. Introduction 

A profound transition in our understanding of stars is now underway, and one of the major 
drivers is the detection of pulsations in large samples of sun-like stars. Much of the immediate 
interest and effort has focused on using scaling relationships for pulsations to infer global properties, 
such as mass, radius, and age. Many new insights into stellar structure will in fact emerge from 
our new ability to design experiments: for example, using masses outside of binary systems, or 
ages outside of star clusters. However, our deepest insights are likely to emerge from diagnostics 
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of internal structure uniquely available from seismic data. In this paper we focus on one such 
property: the depth of the surface convection zone (CZ). We demonstrate that theory predicts 
strong mass and composition trends in the depth of the CZ with surprisingly small errors. This 
raises the prospect of an absolute seismic abundance calibration and rigorous tests of interiors 
theory. Furthermore, we may be able to test more challenging issues, such as deep mixing in the 
mid-F star lithium dip, or the mixture of heavy elements. 

A rich spectrum of non-radial pulsations is observed in the Sun. The study of solar oscillations, 
or helioseismology, has yielded insights into solar and stellar structure. In the case of the Sun 
we have spatially resolved information and can reconstruct the speed of sound as a function of 
depth. For stars we can only observe global modes 0, which makes full sound speed reconstructions 
impractical. The low / modes have asymptotic frequency spacings predicted by theory, however, 
which can yield valuable information about the global properties of stars. Sharp localized changes in 
structure will also manifest themselves as deviations from these regular spacings; examples include 
ionization zones and transitions from radiative to convective energy transport. We briefly review 
the former before introducing the latter, which are the main focus of the current work. 

The average relationships between frequencies of mode pairs {l,n) — {l,n — 1) and (/,n) — 
(Z + 2, n — 1) are frequently referred to as the large and small frequency spacings, respectively. 
The former is related to the mean density, while the latter is a measure of the degree of central 
concentration, and thus he lium content , of main sequence stars. The small frequency spacing 
is a potent age diagnostic (jUlrichI 119861 ). A significant advantage of these relationships is that 
useful information can be extracted from the average differences of many mode pairs, increasing 
the effective signal to noise. The frequency of maximum power reflects a competition between the 
spectrum of turbulence generating the sound waves and the acoustic cutoff frequency. The cutoff 
frequency is empirically observed to follow regular scaling relationships, and the combination of 



Vmax and Au can be used to solve for the mass and radius (Brown et al. 



Belkacem et al. 



2011 



Chaplin et al.ll2011al ). 
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Mosser et al 



201C 



With better data it is possible to extract entirely new kinds of information. It was rec- 
ognized early on that the solar oscillations could be used to measure the stren gth of the he - 
lium ionization zone, and thus by extension the solar surface helium abundance (jCoughl 11984 ) ; 
Delahave &: Pinsonneaultl (|2006l ) found a literature average of 0.248 it 0.004. The surface helium is 
significantly below the initial levels, of order 0.27, required to reproduce the solar luminosity at the 
solar age. This differenc e can be attributed to gravitational settling of he l ium and heavy elein ents 



(|Aller Chapman 
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Noerdlinger 
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The transition from convective to radiative energy transport induces a discontinuity in the 
temperature gradient, which in turn produces a characteristic response in the frequencies of the 



^Modes are specified by the three spherical harmonic quantum numbers: n, the overtone, I, the degree, and m, 
the azimuthal order. For stars other than the Sun, we can generaUy detect only modes of low degree < 3). 



- 3 - 



modes which cros s this boundary. This phenomenon can be used to infer the depth of the surface 
convection zone (jChristensen-Dalsgaard et alj Il99ll ) . The agreement between the theoreticaUy 
predicte d depth and the seismic d ata w as a significant triumph for stellar interiors theory; see for 
example iBahcall &: Pinsonneaultl (|l992l ) . 



In the solar context these seismic tools can be used for precision tests of astrophysics; for 
example, the convec tion zone depth can be measured to a remarkable precision, of order 5 x 10~^ 
(|Basu &: Antial 120041 ') . Scalar constraints on convection zone dept h and surface helium can be corn - 
bined to test the absolute solar metal abundance an d mixture (iDelahaye &: Pinsonneaultl l2006l) 



19891 ) can reduce the effects of settling ([Richard et alJ Il996l : iBahcall et al 



Bve 

Mixing sufficient to explain the solar lithi um depletion (IWeyrnann &: Seardll965l: 



Pinsonneault et al, 



2001 



) and has a char- 
acteristic signature in the sound s peed profile. Solar data is even of sufficient quality to detect 
the signature of metal ionization (iBasu &: Antial 120081 ). which can be used to infer the absolute 
oxygen abundance. In stars i t is anticipated th at the acoustic radius of the convection zone can 
be measured to of order 2 % (iBallot et al.ll2004l'l. M easurements of the surface helium abundance 
are more promising in red giants (jMiglio et al. 2010), but will be technically more difficult in main 
sequence stars; similar comments apply to the extent of convective cores. We therefore focus this 
study on the surface convection zone depth. 

To first order the depth of the surface c onvection zone is deterni ined by the effective tem- 
perature, modified by the surface gravity (see iPinsonneault et al.ll200ll . for a discussion). Stellar 



metallicity plays a second order role, in the sense that lower bulk metallicity implies a shallower 
surface convection zone at fixed effective temperature. Our first priority is therefore to quantify 
these effects, and determine the robustness of the theoretical predictions. Mixing and element sep- 
aration processes can also modify the convection zone depth at a detectable level, although it will 
clearly be more difficult to detect these more subtle shifts in the seismic properties. We therefore 
also explore our ability to empirically diagnose these phenomena. 

We begin with a discussion of our methods in Section [2l We discuss our models results, both 
in terms of the overall trends in the depth of the convection zone, as well as their magnitude in 
comparison to theoretical and observational uncertainties in Section [3l We discuss further potential 
uses of the depth of the CZ as a diagnostic and conclude in Section HI 



2. Calculation of the Acoustic Depth and Theoretical Errorbars 

Our overall approach is similar to that employed in solar model studies. We define a reference 
model calculation which includes our current best estimates of the input physics. We then calibrate 
a solar model and run a series of models with different masses and compositions; these define our 
predicted theoretical trends. We then perform a comprehensive error analysis to estimate the 
detectability of these signals. This includes both theoretical errors (for example, nuclear reaction 
cross-sections or quantum mechanical opacity calculations) and observational errors (for example, 
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the seismic age and mass of the star being tested). Errors are not easy to estimate in some 
cases, such as for opacities. We therefore use the differences between competing calculations as 
our measure of these theoretical uncertainties. We also explore some indirect tests of physics not 
typically included in interiors models, such as rotationally induced mixing. In a full calculation 
one would properly account for nonlinear effects rather than employing a strict parameter variation 
study. We briefly discuss some such cases, but our main concern is with laying out the baseline 
theoretical expectations and their overall reliability. More complex nonlinear calculations would 
be a logical next step and would be better motivated in the presence of a significant observational 
database. 



2.1. Calculation of the Acoustic Depth to the Convection Zone 



We choose to focus on the acoustic depth of the convection zone, Tcz, rather than its physical 
depth, sinc e r^z is the a steroseismic observable of the CZ location. In general, the "acoustic depth" 
is defined (|Goughlll99d ) as, 

fR2 ^f>l 

— , (1) 

Ri Cs 

where R is the radius and Cs is the sound speed. In order to easily compare Tcz for stars across 
a wide range of masses, we use the "normalized" acoustic depth, t^z n, given by. 



j-R* dR!_ 



(2) 



R=0 



where the sound speed is simply Cg = y 

Both P and p as a function of radius can be obtained directly from a converged stellar model. 
In general, the adiabatic exponent, Fi, is a combination of the quantities 



dlnP 
dlnp 



-Xp, 



(3) 



ad 



where Cy and Cp are the specific heats at constant volume and pressure, respectively, and are related 
by 

(4) 



where the derivatives xt and xp are defined as 



XT 



Xp 



pT Xp 



dlnP 



dlnT J p 

dlnP^ 
dlnp J J. 



(5) 
(6) 



- 5 - 



In the case of our stellar models, the values of Ti come either from tabulated equations of state, 
or from a combination of Cp and the derivatives XT and xp: which can be calculated within the 
evolution code itself. 



2.2. Standard Input 



We define a standard set of input physics from which we create our fiducial stellar models. 
Models include both helium and heavy element diffusion, since heavy elements sink with respect 
to lighter elements in a gravitational p otential. We us e the procedure of iThoul et al.l (| 19941 ). which 
numerically solves the full set of the iBurgera (j 19691 ) equations for a multicomponent fluid with 
no restriction on the number of species considered. For the purposes of computin g diffusion co- 



efficie nts, we treat all heavy elements in the same manner as fully ionized iron (see iBahcall et al 



19951 ). In reality, the effects of diffusion and settling are modified by at least two other physical 
processes: radiative levitation and mixing. The first is a small effect for the mass ranges we con- 
sider, and our models therefore include no prescription for levitation (see Section [J] for discussion). 
We do, however, account for mixing. Mild enve lo pe mixing is needed to e xplain the Li and Be 
abundances of low mass stars (jPinsonneaultl 119971 ). iRichard et al.l (| 19961 ) and lBahcall et al.l (120011 ) 
found that mixing sufficient to explain the observed Li depletion has the primary seismic effect of 
re ducing the efficiency of element segregation. We therefore set the diffusion coefficients to 0.8 as 



m 



Delahave &: Pinsonneaultl (|2006l ) to account for the effects of mixing. 



The atmosphere and surface boundary conditions are given by the iKuruczj (|l997l ) model at- 
mosphere t ables H. The convection zone depth is only weakly sensitive to the choice of boundary 



conditions (Bahcall &: Pinsonneau. 



I992I : IPelahaye &: Pinsonneaultl I2OO6I ). We choose the Kurucz 



tables over those of lAllard et al.l (|2000l ) because of their finer gridd ir ig in composition. We uti- 
lize the recentl y updated nuclear reaction rates of lAdelberger et al.l (120111) with weak screeniri g 
( Salpeter 1954 ). and employ the mixing len gth theory of convec tion ( Cox 19681: Viteris e 'l953V 
Opacities are from the Opacity Project (OP) ([Mendoza et al.ll2007l ) for a lGrevesse &: Sauvali (|1998i . 
hereafter GS98) solar mixture, and are interpola ted for each cornpositio n as needed. These are sup- 
plemented with the low temperature opacities of lFerguson et al.l (j2005l ). also fo r the GS98 mixture . 



The GS98 solar abu ndances are in good agreement with asteroseismolo gv fsee iBahcal. 



Basu &: Antiall2008l ) in comparison to the more recent solar mixtur e of lAsplund et al 



et al. 



are thus our default choice. We discuss the effects of adopting the lAsplund et al.l (I2OO9I )" mixture 
in Section 12.3.41 The structural effects of rotation and convective overshoot are neglected in the 
standard models, although we address such processes further in Sections 12.3.21 and HI Semiconvec- 
tion (see iKippenhahn fc Weigertlll994l ) is nominally included, although is of little importance over 
the stellar mass range we consider. 



2005 



1120091), and 



^Models are available at 



http: / /kurucz. harvard.edu/l 
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We utilize the updated 2006 OP AL equation of state (EOS fl (|Rogers et al.lll996l : lRogers &: Navfon"^ 



2002) and the lSaumon et al.l (|l995l ) EOS for temperature and density combinations outside of the 
OPAL tables. To calculate the sound speed throughout each stellar model in a thermodynamically 
consistent fashion, we use the published values of Fi in the OPAL 2006 EOS. Pi is determined 
using the values for P and T from the converged model at the end of each timestep, for each shell 
of the model. Values of Pi for the envelope and atmosphere are likewise acquired directly from 
the OPAL EOS tables, using the values of P and T from the envelope integration. For our pur- 
poses, we neglect the acoustic thickness of the atmosphere in the calculation of Tcz,m because the 
modes of interest are generally evanescent in this region (but see Section [3.51 for further discussion 
of the atmosphere). The base of the convection zone, Rcz, is defined to be the location where 
the Schwarzschild criterion for convective instability, V^ad > V^^, is fullfiled. With the values of 
Fi obtained from the OPAL 2006 EOS, and the structure from the interior and envelope calcula- 
tions within the code, we perform the integral Tcz,n by interpolating the calculated values of Cs[R) 
onto an even grid in ii, and the integrating the tabulated values using a five-point Newton-Cotes 
integration formula. 

We use a solar calibration to set the value of the mixing-length parameter, a (the ratio between 
the convective mixing length and pressu re scale height), and the initial composition, X, Y, and Z 
such that a IMq model at 4.57 Gyr (see lBahcall et al.lll995l ) recovers the solar radius, luminosity, 
and surface abundance of Rq = 6.9598 x 10^° cm, Lq = 3.8418 x 10^3 ergs s'^, and Z/X = 0.02289 
from GS98, respectively. A calibration using this standard set of physics yields a = 1.93271, 
X = 0.710040, and Z = 0.018338. 



2.2.1. Composition grid 

We created a larger grid of models for masses 0.4 Mq - 1.6 Mq and initial abundances —1.2 < 
[Z/X] < +0.6 where [Z/X] = logio ((Z/X)^o,e//(Z/X)0,i) and {l/X)^- = 0.025828, as opposed 
to the GS98 surface abundance of 0.02289 (the difference being due to element diffusion). We 
normalize to this initial solar Z/X throughout the paper, which amounts to a zeropoint offset of 
0.0524 dex between a metallicity scale normalized to initial versus surface solar abundances. We 
use models with the standard set of physics to investigate the effect of composition on the location 
of the convective boundary. The mass range is chosen to roughly coincide with the onset of fully 
convective models on the low-mass end and vanishingly thin convective envelopes on the high-mass 
end. The choice of metallicities is motivated by the typical distribution we expect to observe in a 
sample of field stars. Models are evolved until they leave the main sequence, or until 14 Gyr has 
elapsed, whichever occurs first. For stars with M > 1.3M0, the convective envelope becomes less 
massive than the default fitting point (1.24 x IO^^Mq) between the interior and envelope solutions. 
The fitting point is moved to a minimum mass of 1 x IO^^Mq to accommodate these models. The 



''updated 2006 tables available at [http: / /opalopacity.llnl.gov/opal.html| 
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grid is composed of models spaced every 0.02 Mq in mass, 0.2 dex in [Z/X], with initial helium 
mass fractions of 0.24, 0.26, and 0.28 and Yi^Q for each combination of mass and metallicity. The 
result is a grid of some ~ 2400 models. 



2.3. Theoretical Errorbars on the Acoustic Depth 

We examine the theoretical errorbars on the acoustic depth through comparisons of pairs of 
model grids. The first grid contains "standard" models in the sense that they represent the results 
for the set of input physics described in Section [2.21 Comparison grids are identical to the standard 
grid except for a single alteration to the input physics. Both grids are subject to separate solar 
calibrations. We divide the parameter variations into several distinct classes, based on the nature 
of variation. Some parameters, such as the diffusion coefficient or nuclear cross sections have well 
defined and random errors. Changes to other parameters, such as the EOS and opacities, can 
also shift the location of the convection zone in either direction, but the uncertainties in Tcz,n 
incurred from switching between different EOS or opacity tables are systematic in nature. We 
treat the changes induced i n Tr.y.^n due to well-motivated y ariations of these parameters as effective 



2a errorbars on Tcz,n (as in lBahcall &: Pinsonneaultlll992l ). There are also uncertainties that arise 



simply from our inability to measure the mass, radius, composition, and age of real stars with 
perfect accuracy, which we will describe as "observational" in nature. A final class, which we will 
term "zeropoint uncertainties" is related to assumptions such as the heavy element mixture or 
the presence of rotational mixing, and for these cases, the resultant theoretical uncertainties are 
asymmetric. For example, one can make Tcz,n smaller by including mixing, but never larger. In these 
cases, we incur uncertainties in the zeropoint of our relations between various physical parameters 
and Tcz- These are considered separately from the systematic, random, and observational error 
sources throughout the paper, but are discussed here for completeness. We proceed, then, to 
address the uncertainties from each of these sources in turn. The values of X, Y, Z, and a for each 
of the physics variations are listed in Table [H 



2.3.1. Random uncertainties 



Element diffusion in stars allows heavier atoms to sink relative to lighter ones. The effect of 
diffusion is to situate metals, which are a significant source of opacity, deeper within the star than 
they would otherwise be, resulting in a deeper convective boundary than in mo dels with no element 
diffusion. The presence of element diffusion in the Sun produces a 1.7% effect feahcal l et al 



on th e location of the convection zone in Solar models (see also iBasu et al.l .2000, : .Bahcall et al 



2001 



20041 ). We const ruct a calibrated s et of models with the helium and heavy metal diffusion coefficients 
altered by 15% (jThoul et al.lll994l ). to mimic uncertainty in the strength of diffusion in the interior. 
Apart from the differing solar calibrations and adjustemnt of the diffusion coefficients, these models 
are identical to those run with standard physics. 
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We address the effect of the nuclear reaction rates on the stellar stru cture, and adopt error 
estimates for nuclear reaction cross-sections from lAdelberger et al.l (|201ll ). The major reactions 
considered are the primary pp chain reactions Si^i (pp), 83^3 (He^ + He^), 83^4 (He'^ + He^), and 
CNO cycle 81^14 (P + N14), which were each changed by ±4(7. 



2.3.2. Systematic uncertainties 



We chose a different prescription for the equation of state in an effort to quanti fy the change 
in Tcz,n due to quantum mechan ical uncertainties. We us e the l8aumon et al.l (jl995l ) (8CV) E08 
instead of the OPAL 2006 E08 ([Rogers &: Navfonovll2002l ) chosen for our standard set of models. 
We choose this particular variation since the differences between earlier versions of the OPAL EOS 



are s mall (see iBahcall et al.l 12004 ). and the relative simplicity of the Yale E08 (jGuenther et al 



I992I ). which treats the interior as fully ionized and solves the 8aha equation for the envelope, 
is a poor representation of the the state of the art to which E08 calculations have progressed. 
The OPAL and 8CV equations of state are both sufficiently modern, and yet have very different 
approaches to the problem, and so the 8CV E08 serves as a useful comparison. Because we draw 
the value of Pi in the standard models directly from the E08 tables, we must alter the manner in 
which we calculate Pi for the 8CV E08 models. Using the relationships between the derivatives 
Xpi XT and the specific heat at constant pressure, Cp (Equations [3] HI EJ [6]) , calculated numerically 
within YREC, we combine these values to calculate Pi and the sound speed throughout the star. 

There ar e two primary sour ces for high temperature opaci ties in stellar interio rs models, the 
OPAL group (IRogers et al.lll996l ) and the Opacity Project (OP) (jBadnell et al.ll2005l ). each of which 
approaches the quantum mechanical calculation of the high-temperature opacities in a fundamen- 
tally different way. A thourough comparison , and discussion of the differences between the two 
methods is present in l8eaton &: Badnelll (|200J). The differences in the Rosselan d mean opacity for 
the c onditions found at the base of the solar convection zone are of order 5% (|8eaton &: Badnell 
2004 ). The opac ity plays a significa nt role in determining the precise location of the base of the 
convection zone ([Bahcall et al.ll200ll ). and so we test the sensitivity of Tcz,n to our choice of opacity 
table by running variant models using the OPAL opacities, instead of our default choice of the OP 
opacities. 

We expect that the choice of atmosphere and boundary conditions will be most important for 
very cool stars. We test the worst-case dependence of Tcz,n on the choice of atmosphere boundary 
condition by creating a calibrated grid of models for a grey atmosphere boundary condition. We 
note that this exercise only quantifies the dependence of Tcz,n on the boundary condition, since the 
portion of Tcz due to the atmosphere, Tatm-, is neglected in the calculation of Tcz,n- 



We also consider the importance of convective core overshoot (see IZahnlll99ll : lMaedeiill975l ) 
to the determination of Tcz.n- While in principle overshoot in all convective layers is possible, 
and affects the local composition, we consider convective core overshoot in particular, because the 
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added fuel supplied to the core through overshoot-related mixing could have broader impacts on 
the physical structure of the star. The addition of core overshoot should primarily affect the high 
end of our mass range, where models begin to develop convective cores. We choose a core overshoot 
parametrized by the pressure scale height, with a value of 0.2 pressure scale heights. We do not 
enforce an adiabatic gradient in the overshoot zone or add envelope 'undershooting'; overshooting is 
treated purely as 'overmixing', and the thermal structure is left unchanged. Observationally, were 
significant overshooting to be present, observers would see it as an effective change in the location 
of the convective boundary. Depending on the nature of the overshooting, this could manifest itself 
as either a zeropoint or mass-dependent shift in the location of the CZ. Here we consider only the 
manner in which overshoot induced mixing affects the evolution because of increased fuel supply 
to the core. 



2.3.3. Observational Uncertainties 



We can expect uncertainties in stellar parameters such as T^fj, M, R, p, Tcz,ni Y , and the age 
simply from the nature of the observations with which they are determined. We likewise consider 
our assumption of a particular mass-radius relationship as a form of observational uncertainty. 
In all cases, these are uncertainties dictated by our ability to measure stellar properties, and are 
therefore of a fundamentally different nature than the uncertainties described above. We also note 
that these uncertainties are subject to potential rapid improvement, typically on a time scale shorter 
than those for improvements in opacity tables or equations of state. 

Asteroseismic age diagnostics are sensitive to the helium fraction in the stellar core, a nd thus 
provi d e information a bout how far along its main sequence lifetime a star has progressed (lUlrich 
19861 ). ICreeveyl (j2009l ) suggest s that we can obtain the ages of main sequence stars to within 10% 
of the MS lifetime, and already iMetcalfe et al.l (120101 ) present an asteroseimic age for KIC 11026764 
accurate to 15% with currently existing Kepler data. Future missions, such as GAIA, which aim to 
attain precise parallaxes on a large sample of stars, may eventually allow us to better constrain the 
age based on an absolute luminosity, but we proceed with the assumption that age can be measured 
to this 10% accuracy, and propogate these uncertainites through our models. 

The abundance of helium in stars is notoriously difficult to measure directly and represents 
another source of observational uncertainty. Since luminosity is also a function of the helium 
content, constraints on the mass and the lumi nosity (at fixed and age) lead to constraints 



on the helium, an idea that goes as far back as ISchwarzschildl ( 19461 ). If we take L = AirR'^aT^jj 



and assume reasonable measurement uncertainties in R and T^fj, the uncertainty in L is (t| 



2aR 
R 



+ 



and the uncertainty in Y due to that in L is Cy = (t| 



dY\' 
'dLj 



. Finally, the 



uncertainty propagated to Tcz,n is then cj^ = cr| 



dr 



'Y , QY^ 

from our composition grid. We assume that R can be measured to a fractional uncertainty of 2% 



. We calculate these derivatives numerically 
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and T^ff to within lOOK. 

A portion of the observational uncertainties wih come directly from the asteroseismic mea- 
surements themselves, namely our ability to precisely constrain the large frequency separation and 
glitch signatures. We assume that t he value of p is atta inable to within a relative uncertainty of 1% 
from asteroseisrnic me asurements (IVerner et al.ll201ll ). and that Tcz,n can be measured to within 
2% (|Ballot et al.ll2004h . 

We have assumed for our standard grid that a single mixing-length parameter, calibrated for 
a solar model, is valid over the entire range of masses and compositions we consider. Because we 
employ only the mixing length theory of convection in these models, we have no ability to test 
how Tcz,n changes as a result of different theoretical assumptions about convection, and rather we 
choose to view variations in a as an uncertainty in the mass-radius relationship (but see Section H] 
for a further discussion of convection theory). Because the general approach in stellar modeling is 
to determine the mass and calibrate the model such that the correct radius is recovered, we treat 
this as an observational error. To test how severely this may impact the inferred depth of the 
CZ, we construct a grid of many different values of a, and choose a(M) such that for all masses 
considered, the radius of the star in the "altered physics" grid is about 1% larger than in the single, 
solar calibrated a case in the standard grid. Observed discrepancies from the theoretical mass- 
radius relationship are observed to be as large as 1 0% for cool , low- i nass stars in binary sy stems 



where they can be well studied (IKraus et al 



2011 



Irwin et al.ll2009l : iBavless fc Oroszl l2006l ;i . On 



both the high and low mass ends of the mass range, the stellar radii become rather insensitive to 
changes in the mixing-length parameter. For high mass stars, this is because the pressure scale 
height is small enough that large changes in a itself are physically of little significance. On the low 
mass end, because stars are nearly fully convective, changes in a tend to shift stars along the main 
sequence, rather than changing the relation. In both extremes, no change in a is ever sufficient 
to produce a 10% difference in the radius. Our chosen AR = 0.01 is achievable over nearly the 
entire mass range considered with sufficient changes to the value of a. Although scaling these 
theoretical errorbars to larger uncertainties in the stellar mass-radius relation is not unreasonable, 
it is important to note that for stars with M < 0.6Mq and M > IAMq, the model radius becomes 
insenstive to a and simple scalings will fail. 



2.3.4- Zeropoint uncertainties 



Observational and theoretical evidence suggests that some form of mixing operates in both 
the Sun and other stars, and that this mixing can have effects on the apparent efficiency of dif- 
fusion. We know from modelling of the Sun that diffusion alone does not adequately repr oduce 
solar light-element depletion relative to meteorites ( Richard et al. 1996 : Bahcall et al.|[200ll ). and 
that rotationally induced mix ing provides a well-motiv at ed physical process by which the observed 
depletion could be achieved (jPinsonneault et al.lll989l ). iBalachandranI ()1995l ) finds that diffusion 
alone cannot explain the Li abundances in M67, and it is generally believed to be the signature 
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of some form of deep mixing. iDelivannis et alj (|l998l ) likewise finds correl ated Li and Be deple- 
tion patterns in Hyades F-stars (the strong Li depletion first recognized bv lBoesgaard &: Tripicco 
19861 ) is best supported by a slow mixing of stellar material. The primary consequence of mixing 
for our purposes is the accompanying decrease in the efficiency of element diffusion, possibly to 
the p oint to which it a ppears that diffusion does not op erate at all. Recent observations of NGC 



6397 (IKorn et al.ll2006l ) and theoretical modelling efforts (IChaboyer et al. 



1992 



Potter et al.ll2nn8l ) 



found that diffusion must be partially suppressed in order to explain the observed trends in the 
light elements of metal-poor stars. To mimic the effects of strong mixing, we construct models with 
no helium or heavy element diffusion. 

The mixture of heavy elements, even from the Sun, is another important systematic error 
source. We consider here both the case of revised Solar heavy element abundances, and the case 
of a-element enhancement in low-metallicity models. In both cases, the relative abundances of 
important contributors to the opacity are altered with respect to iron, and we seek to investigate 
the sensitivity of Tcz^n to changes in the element mixtures.. 

In the case of the Sun, there exists a well known tension between the solar CZ depth and sound 
speed profile inferred frona helioseismology versus that implied by the recent solar abundances of 
Asplund et al.l (j2004l . l2009l ). These recent abundances are based on non-LTE, 3D radiative trans- 
fer calculations of the solar atmosphere and represent the state-of-the-art in modern atmosphere 
modeling. However, solar models constructed with this new, low bulk metallicity are in worse agree- 
ment with seismic diagnostics such as the surface heliu ni abundance, CZ d e pth, and solar soun d 



speed profile than models with the older GS98 mixture (jBahcall et al 



2005 



Basu fc Antiall2008l ). 



Although the new solar mixture has a similar iron abundance as the old GS98 mixture, the CNO 
elements are significantly adjusted, and the oxygen abundance in particular is quite low. Because 
these elements tend to be completely ionized in the deep interior of stars, they contribute little in 
the way of opacit y in the core, but have significaii t opacities near the location of the base of the CZ 
in solar-like stars foelah ave fc Pinsonneault 2006 ). It is important to note that similar work on the 



Solar mixture bv.Caffau et al 
metallicity than lAsplund et al 



201 ih. also using a 3D analysis, arrived at a higher oxygen and bulk 
(120091). Although fut ure work may alleviate the conflict between 



helioseismic inversions and the lAsplund et al.l (120091 ) mixture, we choose to investigate how the 
lowest published oxygen abundances affect our conclusi ons regarding the dep th of the convection 
zone. We con stru ct a calibrated s et of models for the lAsplund et al.l (j2009l ) mixture, using OP 
(jSeatonll2005l ) and iFerguson et al.l (j2005l ) low temperature opacity tables adjusted for the change 
in mixture. The models are calibrated to a surface abundance of Z/X = 0.0199. We note that 
while the pre-main sequence stellar models and opacity tables are adjusted to reflect the difference 
in abundance pattern, the EOS and atmosphere tables are not. In both cases the correct mass 
fraction in metals is used, and the errors incurred from the difference in mixture in the atmosphere 
and EOS should be negligible, since the change of mixture primarily affects the CNO elements and 
therefore nuclear burning and the highly metal-sensitive opacities. 



In the case of metal-poor halo stars, we may also expect that there may be deviations from 
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the Solar mixture due to differ ent chemical enrich ment histories. We consider a-enhanced models, 
with [a/Fe] = +0.2 (following iDotter et al.l 120081 ). We compare models with the standard GS98 
mixture at [Z/X], [Fe/H] = -1.0 to a-enhanced models with [Fe/H] = -1.0 but [Z/X] = -0.85 
with Yi = 0.247 in both cases. As with the Solar case, both low and high temperature opacity 
tables and initial models are adjusted for the change in mixture. We supply EOS and atmospheres 
with the correct bulk metal abundances, but relative abundances are not adjusted. 



2.3.5. Combined Uncertainties 

We calculate the magnitude of the uncertainties from each of the error sources above for an 
assumed reference model at solar composition, with an age of 5 Gyr, with the standard set of 
physics. In combining the errors from each source, we treat the uncertainties as uncorrelated. We 
add individual sources of error in quadrature. For example, the random error on Tcz,n-, crr,rand is 
composed of cr'^j.^nd ~ ^'^nuc + '^diff^ where the individual errors are from the nuclear reaction 
rates and diffusion coefficients, respectively. Similar combinations of error terms are calculated for 
each class of uncertainty, random, systematic, and observational. The total uncertainty on Tcz,n is 

2 _ 2 _|_ 2 I 2 

where the random, systematic, and observational errorbars are added in quadrature. We consider 
the zeropoint uncertainties separately, since they are of a fundamentally different nature, and have 
asymmetric effects on the depth of the convection zone. 

Although an exhaustive investigation of the cross terms in our theoretical error estimates are 
beyond the scope of this paper, we do comment briefly on a few cases in which we have investigated 
the role of crossterms with the age uncertainties. The shapes of the curves in Tcz,n—Teff and Tcz,n—p 
space depend strongly on age, and so we have singled out this particular error source in which to 
look for crossterms. We find that for the case in which both the age and helium are simultaneously 
considered, that the error on Tq^^^ from the combination of age and Y uncertainties is negligible 
when the interplay between the error sources is considered. Crossterms between age uncertainties 
and changes in the physics are likewise negligible. The exception is in the case of overshoot models, 
in which the time dependence of the diffusion and the age uncertainties interact, inflating the 
errorbars by up to 30% for the high mass models when both age and overshoot uncertainties are 
considered simultaneously. This is not unreasonable, since overshooting is both more important for 
more massive objects, and affects the amount of fuel in the stellar core. In general, we recommend 
that one estimate errors, and potential cross-terms, on a source- by-source basis when actual data 
are available. 
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3. The sensitivity of the convection zone to mass and composition 

In this section we present the main predictions of our models, the results of which are given 
in Table [2j The depth of the convection zone is a strong function of mass and composition, as 
predicted by interiors theory and recovered in our models. With our theoretical errorbars, we find 
that these trends are a strong, observable effect, even when we account for uncertainties in both the 
interiors models and observationally derived quantities. The mass and composition dependencies 
can be viewed in terms of the purely asteroseismic quantities of the mean density {p, here given 
as p = M/B?) and normalized acoustic depth Tczl^i, = Tcz,n, and any analysis can in principle 
rely on solely asteroseismically obtained measurements. Furthermore, the composition dependence 
of Tcz,n is more pronounced in p - Tcz,n space, as opposed to, for example, T^fj - Tcz,n space: the 
use of purely asteroseismic variables is not only useful, but beneficial. Finally, we show that with 
appropriate observations, various tests of stellar physics, beyond the basic question of the location 
and presence of a convection zone, are possible. 



3.1. Mass dependence of r, 



cz,n 



We expect from very simple interiors arguments that the depth of the convection zone must 
be a strong function of mass. To first order the location of the base of the convection zone is 
set by the location of the H and He ionization zones, where the adiabatic temperature gradient is 
suppressed below the radiative temperature gradient and the criterion for convection is satisfied. 
This is illustrated in Figuredl which shows Rcz/Ri, versus mass on the left, and Tcz,n versus p in the 
middle and Tcz,n versus T(.fj on the left for standard models at solar composition and an age of 1,5 
and 10 Gyr. The expected strong mass dependence of the depth of the CZ is clearly present, and 
the mapping from the Rcz/R* vs. M to Tcz,n vs. p planes is a simple one, which preserves the sense 
of the trends. Models with relatively high masses have vanishingly thin, shallow convection zones, 
whereas low mass models are nearly fully convective. If our physical understanding of what sets 
the location of the convection zones is corre ct, observations of rn any different stars should display 
this strong predicted trend (also discussed in lMonteirn e,, .1 I B In the case in which this strong 
dependence is not observed, we immediately learn that there is some fundamental physical process 
that has been neglected or poorly treated in interiors models. 

We can also comment here on the basic time dependence of such relationships. At very young 
ages, the entire mass range we consider (0.4 < Mq < 1.4) is on the main sequence. The basic 
shape of the Tcz,n vs. p changes very little between 0.5 and 1.0 Gyr, for example, because the 
most massive objects we consider have main sequence lifetimes of a least 1 Gyr. Once we begin to 
look at later times, however, the low density tail of the curves begins to show significant changes, 
because the stars that occupy that part of parameter space are progressively less massive (more 
massive models have evolved off the main sequence and out of our realm of consideration). The 
mean density of a star decreases over the course of its main sequence lifetime, so it is possible for 
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an older, less massive star to have the same mean density as a young, more massive star. Their 
convection zones, however, generally not at same relative depth on the main sequence, and so the 
low density tail of a given Tcz,n vs. p curve shifts to deeper convection zones at late times. We 
incorporate this feature of the relation into our observational error estimates in later sections. 



3.2. Composition dependence of n 

While mass should be the primary determinant of the depth of the convection zone, it is clear 
from similarly simple arguments that the composition should also p lay some role in the location o f 



the convection zone. The radiative temperature gradient is given bv lKippenhahn &: WeigertI (|1994| ) 
as Vrad = i67raeG where o, c and G are the usual physical constants, I is the luminosity of the 
shell with mass m, P and T the pressure and temperature, and k the opacity. With all else being 
equal, an increase in the opacity leads to an increase in the radiative temperature gradient, which 
in turn means that the criterion for convection is satisfied deeper (at higher T) within the star. 
While metals are an almost negligible fraction of the mass, because they contribute significantly 
to the opacity and are important in determining the precise location of the convection zone. This 
is indeed what we find, shown in Figures I2|31 HI and [5l where Tcz,n and Tcz (in seconds) is plotted 
with respect to p and Tg//. We predict that Tcz,n changes by 0.5-1% per 0.1 dex in [Z/X] over most 
of the mass range we consider. Furthermore, this composition signature is an absolute, rather than 
relative, measure of stellar abundances. It is most pronounced in the Tcz - p plane, rather than in 
the Tcz ' "^eff plane. Therefore, the most preferable space in which to work is also the space in 
which the composition measurement can be made solely with asteroseismically obtainable variables: 
spectroscopic and photometric characterization of the stellar parameters is only necessary as an 
additional constraint on the stellar parameters. 

Because our models include gravitational settling, heavy elements tend to sink relative to light 
ones, and the surface [Z/X] is generally not the same as the initial abundance, and changes as a 
function of time. Figure [6] shows both the difference between the initial and surface abundances. 
This difference arises because of gravitational settling of heavy elements and would manifest itself 
as a Tcff dependence of the surface [Z/H] in a sample with homogeneous initial composition, such 
as an open star cluster. If we consider models at fixed p/pq = —0.2 at 1.0, 5.0, and 10.0 Gyr as in 
the right panel of Figure El we find that the difference between the surface and initial abundance 
is most pronounced for the 10.0 Gyr curve. This is a balance between two competing factors: 
at earlier times, more massive stars with short settling timescales are still on the main sequence 
and occupy this density range. At later times, less massive stars occupy this region and have 
longer settling timescales, but longer MS lifetimes over which settling can occur. This difference 
between the initial and surface abundances is important for any comparison of asteroseismic and 
atmospheric abundance measurements: the value of the surface abundance for a given model is a 
physics and age dependent property. 

We also show in Figure [7] the fractional difference in Tcz n among models of different initial 
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helium abundances at constant, solar Z/X at 1.0 Gyr. We factor the uncertainty in the helium, as 
discussed in Section [2.3.31 into our observational error budget. 

3.3. Uncertainties in the relationship of Tcz,nip,Teff) 

The characterization of the uncertainties in the relationships among mass, composition, and 
Tcz,n means that we can not only comment on the existence of an important trend, but also quantify 
whether it is presently observable. We first discuss the results of the variant models introduced 
in Section 12.31 and then show that these uncertainties are small enough that the depth of the CZ 
can be used as a precise indicator of composition. We present representative uncertainties in Tcz^n 
due to random, systematic, and observational uncertainties in Figure [8] for models with a solar 
composition and age of 5 Gyr. 

The contribution to o",- from sources such as the diffusion coefficients and nuclear reaction rates 
is below 0.1%, except on the very low and very high mass ends of the distribution. Contributions 
from the systematic class of errors are somewhat more significant, with the EOS being the most 
important source of uncertainty in this particular grouping. As expected, the uncertainties due to 
overshoot and the choice of boundary condition are small, below 0.05%. 

Uncertainties incurred from errors in our knowledge of the global properties of the star (M, T^ff, R, Y 
and age) are by far the largest source of error in Tcz,n- At low mean densities, uncertainty in Tcz^n 
due to uncertainty in the age is the most significant contributor to the observational error, which 
is unsurprising given the behavior of the curves in Figure [T] as a function of age. At small masses 
the uncertainty in the age itself is large because the objects have very long MS lifetimes and we 
can only measure the age asteroseismically to within 10% of that lifetime. However, the change in 
Tcz,n with time is comparatively small and so the uncertainties in Tcz,n induced by age variations 
are modest for the low-mass models. For higher masses, our ability to measure the age asteroseis- 
mically is substantially better, but the shape of the the Tcz,n vs. p is changing significantly with 
time because models are evolving off the main sequence and a given density probes very different 
masses at different times. Therefore, the age-induced uncertainties are largest on the high mass 
(low density, high Tg//) parts of the curve. 

One should note that the uncertainties due to the mass-radius relationship are relatively small, 
but that we have also chosen a very modest AR = 1%. In principle, these errors can be scaled 
for larger radius uncertainties in the mid-mass range. On both the large and small mass extremes, 
however, simple scalings of radius errors using the mixing length a will fail (as mentioned in l2.3.3p . 
On the low mass end, even a very large change in a has only a small effect on Tcz,n- For high masses, 
however, large changes in a also lead to substantial changes in Tcz,n- Therefore, one must be wary 
if attempting to scale these particular errorbars for larger radius discrepancies for the higher mass 
stars. 

The uncertainties due to unknown helium abundances are also non-negligible contributors to 
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the observational error. On average, given the assumed uncertainties in our abihty to measure lumi- 
nosity, we can hope to constrain the helium mass fraction to within 0.01-0.02, using the technique 
described in Section [2.3. 3[ Additional information from parallaxes or asteroseismic determinations 
of the helium could better constrain this number. It is however, encouraging that we will be able 
to assign realistic errorbars to the helium, as opposed to ad hoc estimates. 

We combine this suite of systematic, random and observational uncertainties on Tcz,n and 
translate this uncertainty into a measure of our ability to measure [Z/X] through o'^i^/X] ~ 

(d[Z /X] \ ^ 
— j , where the components of (Jtcz,„ are shown in Equation[71 Figure [9] shows C[2/X] 

for each error source, systematic, random, and observational, as well as the combined total error. 
The result is that Tcz,n is actually strikingly sensitive to composition, even when we include rea- 
sonable theoretical and observational errors. The uncertainty in [Z/X] is in the range of 0.15-0.3 
dex over the mass range we consider. 

In addition, we alter the compositions and perform the same comparisons for sets of metal-poor 
and metal-rich models (rather than errors derived for the solar case presented in Figure [9|) . We 
assume a simple chemical evolution scheme of the form 

F = r, + fz (8) 

with Yp = 0.246 and dY/dZ = 1.0. We take models with initial Z/X ratios at a tenth solar and 2.5 
times solar, with the change in Y determined by our chemical evolution assumptions. Examining 
the fractional differences in Tcz,n for sets of models at different initial abundances allow us to check 
whether our ability to determine [Z/X] depends strongly on the composition. We find that the 
situation is quite favorable, with f'"[2/X] ~ ^'^ ~ ^'"^ brightest [Z/X]. = — 1.0 models at 10 

Gyr (representative of a halo star population) when all assumptions about the observational errors 
are identical to those in the solar example. In this case of a metal rich object with [Z/X] - = -1-0.4, 
^[Z/X] ^® similar to that in the solar case. This suggests that Tcz,n remains a good indicator of 
composition across the entire regime of compositions we have considered, provided our assumed 
observational errors remain representative. 



3.4. Probing the physics of stellar interiors with Tcz,n diagnostics 

Our analysis suggests several interesting tests of the conditions that prevail in stellar interiors 
using measurements of Tcz,m beyond the potential to constrain composition and confirm basic 
theoretical predictions of interiors models. Because the depth of the convection zone has some 
sensitivity to the particular physical assumptions of interiors models, we can invert the question 
confronted above and ask: if we can measure Tq^^ti^ and if we can trust our stellar parameters 
derived by means other than asteroseismology (photospheric metallicities in particular), can we use 
the patterns we observe in Tcz,n to infer something about the physics of the interior? In this section 
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we first describe how one could use measurements of Tcz,n as a test of the physics responsible for 
the observed mid-F star Li depletion. Secondly, we outline the manner in which one could use a 
large sample of Tcz,n measurements to constrain the stellar abundance pattern. 



3.4.1. The Li dip 

We observe that Li undergoes a severe de pletion event in stars of roughly 6200 — 6350K 



(jBoesgaard Tripiccd 119861 : iBalachandranI Il995l ) . While the existence of the Li dip is well estab- 



lished observationally, theory has yet to converge upon a mechanism responsible for the effect. Many 
different mechanisms have been proposed to explain Li depletion in both the Sun an d other stars: 



throu gh mixing by waves, mass loss, diffusion, and rotationally induced mixing (see iPinsonneault 



19971 . and references therein, for a thorough discussion). We focus here on rotational mixing and 
imagine the following scenario: if stars undergo an episode of strong rotational mixing at about 
~ Q350K, and that mixing has the effect of completely erasing the element segregation induced by 
gravitational settling and diffusion, then we expect a jump in Tcz^n as the model crosses the Li dip 
boundary and the underlying physical assumptions change. We find that the no-diffusion models 
have values of Tcz^n that are ~ 6% lower than those in the standard model over the temperature 
range of 6200-6350 K for models at 1 Gyr (at later times, few models populate this temperature 
range). We imagine a scenario in which stars are well represented by the standard model curve up 
until the edge of the Li dip, at which point they undergo a strong mixing event which erases the 
effects of diffusion, and the star abruptly jumps onto a no-diffusion model curve. If we consider 
the case in which we observe pairs of stars, one of the Li "peak" at 6200 K, the other in the Li 
"dip" at 6350 K, the important quantity to consider is the slope, ATcz,n/ '^T^.j f of the standard 
models over this tempertaure range, compared to that between a standard model on the low tem- 
perature side of the Li dip, and a no-diffusion model at the high temperature side. We consider 
the slope Atcz^u/ ^T^jf of the standard model and the scatter in that slope present when we in- 
troduce our aforementioned changes to the physics, compared to the slope between the standard 
and no-diffusion models over the dip. In the ideal case in which we have perfect measurements 
of Teff and Tcz^n, and the only uncertainties are theoretical (not observational), then the jump in 
the value of the Li dip is an 80" event. However, when observational errorbars due to 

age, Y, and astereoseismic measurement uncertainties are included on the standard-to-no-diffusion 
model slope, the jump in Tcz^n is significant at the 0.8a level per pair of stars In Tcz,n — P space 
the significance is slightly decreased, due to the fact that the mapping between Tf,ff,p and mass 
changes slightly between the standard and no diffusion cases, and conspires in this plane to make 
the jump less visible. We conclude then, that with a sample of ~ 15 pairs of stars, if a mixing event 
is responsible both for removing the signatures of diffusion and providing the means to deplete Li, 
then a trend in the observed values of Tcz,n should be visible at the 3(T level. 

There clearly exist some caveats to this prediction, the most important of which is that it is 
unclear whether rapidly rotating stars of the sort on the hot side of the Li dip will actually display 
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solar- like oscillations and produce reliable measurements of Tcz,n- The detection and interpretation 
of solar oscillations in rapidly rotating stars is still among the principal challenges in asteroseis- 
mology (|Reesdl20ld ). Furthermore, our models have negelected the structural effects of rotation, 
and we have chosen to use no-diffusion models as a simple representation of rotationally mixed 
stars. Rotationally induced changes to the structure, especially in hot, rapidly rotating Li dip 
stars, may be important. Nevertheless, we provide a useful test of the mixing in stellar interiors in 
the temperature range of the Li dip. 



3.4-2. Variant Elemental Mixtures 



We briefly also consider a-enhanced models, with [a/Fe] = +0.2 (following iDotter et al.ll2008l ) 
at fixed [Fe/H] at 10 Gyr, to mimic an old halo star population. If one were to consider a pair of 
stars at the same p, and assume that the theoretical uncertainties on a-enhanced stellar models are 
the same as in the case of a solar mixture, then their measured values of Tcz,n should be different 
by of order 1.0-1.3cr for p < Pq with the inclusion of observational uncertainties. In general, the 
difference between models with solar and a-enhanced mixtures is Atc^^^ ~ 0.03 — 0.06 for the stars 
with p < Pq aX 10 Gyr, which will also be the most likely objects to be detected in missions such 
as Kepler. Therefore, if one can measure the [Fe/H] values for several pair of stars at the same age 
and mean densities, then it may be possible to distinguish between solar and a-enhanced mixtures 
on the basis of the normalized acoustic depth to the CZ. This is another example of the power of 
pairwise comparison, since the theoretical errors effectively cancel for two stars of the same mean 
density and age, and it is only observational errors that effect the significance of the difference in 



As we discussed in Section 12. 3. 4^ the recent lAsplund et al.l (120091 ) oxygen abundances are 
contentious in part because the revision implies a solar CZ depth that does not agree with as- 
teroseismic measurements. We investigate here whether we can utilize ensemble measurements of 
Tcz,n to learn about the oxygen abundance relative to the total metal abundance of other stars. 
This question is well-posed in an open cluster situation, in which the stars are of uniform age 
and composition, and the stellar parameters are somewhat better constrained than in the case of 
a random field star. The typical difference between standard solar models and Asplund mixture 
models is /\Tcz,n ~ 0.005 — 0.01 for models at solar composition at 1 Gyr. We will focus on a 
sample of stars, randomly drawn from a uniform distribution in 5500 < Tf^ff < 6500, each with an 
average combined observational and theoretical uncertainty of a-r ~ 0.03 (including a 0.3 dex [Z/X] 
uncertainty). The quantity Tstandard - Tobserved, where Tstandard is the acoustic depth that would 
be measured for standard physics, and Totserved is the value measured for stars with an Asplund 
mixture, quantifies the zeropoint offset induced by the difference in mixture. Given our standard 
assumptions about the theoretical and observational uncertainties, careful measurements of Tcz,n for 
a 25 star sample could detect a mixture difference at 3a. One should note that this is an idealized 
example: we've assumed that all stars are exactly the same age and composition, with exactly the 
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same oxygen abundance, and exactly the same physical processes operating within them. In reality, 
we could imagine that some effect, such as mixing, might operate differently in stars of different 
masses, which would dilute the abundance pattern signal, or make it appear anomalously strong, 
depending on the sense of the mass dependence. Both mixing and low relative oxygen abundances 
tend to make the CZ more shallow, and so disentangling the zeropoint offset due to a different rel- 
ative oxygen abundance may in practice be quite challenging. This analysis also relies on a correct 
theoretical zeropoint calibration: our standard physics models must be anchored correctly, because 
any theoretical zeropoint offset could be mistaken as a real signal. It is currently unclear how well 
we can achieve this, as even numerical sources of error become important at the Atc^ „ = 0.005 
level. Nevertheless, this is another potentially useful application of careful measurements of Tcz,n- 
In the solar case, detailed information about the mixture could be obtained through a simultaneous 
measurement of the surface helium and Rcz- We anticipate that a similar approach, if practical, 
will be required in the stellar context. 



3.5. Caveats 

We discuss the caveats to our findings in regards to our treatment of the atmosphere and 
convection theory, and other physically important elements of stellar interiors, such as rotation, 
magnetic fields, and radiative levitation. 

We have neglected the contribution to the acoustic depth from the atmosphere throughout our 

discussion: we have only considered the acoustic depth due to the interior and envelope portions 

of our models. Although all models are run with a Kurucz atmosphere boundary condition (P at 

T = Teff), we perform the calculation of the "acoustic thickness" of the atmosphere using a grey 

atmosphere, to allow us to calculate the necessary integrals as a function of radius. While we have 

already demonstrated that the choice of boundary condition produces small (~ 0.5%) changes in 

the normalized acoustic depth, we expect that there are somewhat larger uncertainties associated 

with the atmosphere itself. We integrate the sound speed in the atmosphere using the assumption 

dr 



that the change in radius is given by dr = (jCoxlll968l ). where r is the optical depth, and k the 

—Kp 

opacity. For models with a grey atmosphere boundary condition, we find that the acoustic thickness 
of the atmosphere, Tatm is 50 > Tatm ^ 250 seconds, with Tatm increasing with increasing stellar 
mass. Tatm is typically 5-7% of Tcz for all but the most massive stars with the thinnest convective 
envelopes, where it is a more significant fraction of Tcz- Tatm is ~ 4% of the total acoustic travel 
time in the interior + envelope regions for all masses. The inclusion of Tatm in the normalized 
acoustic depth can change Tcz,n by up to 20% for massive objects with thin convective envelopes, 
but is typically 5% for stars with < I.OMq A calibrated, standard physics, solar model at 4.57 
Gyr produces a Tcz = 2100s, whe reas the solar value for Tcz, which includes surface and atmospheric 



contributions is ~ 2200 — 2300s (jVerner et al.ll2004l ). which suggests our model results are in good 



agreement with reality. From these arguments, we can reasonably expect that the neglect of the 
atmosphere may result in a few hundred second offset between our models and reality. 
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One can also expect that given the nature of our asteroseismic observables, the derived values 
should not suffer substantially from uncertain surface term corrections. In both the case of the 
mean density (derived from the large frequency separation) and the acoustic depth to the CZ 
(derived from an oscillatory signal in what would otherwise be uniformly spaced frequencies) it is 
only the relative difference in the surface term among the modes that will bias the measurements. 
Furthermore, in models where surface convection re sponsible for the asterose ismic surface terms is 
treated more carefully than in our MLT approach ( Stein Nordlundl 2000 ). the inclusion of the 



additional physics significantly improves agreement with high frequency solar modes, but the solar 
base of the CZ remained essentially unchanged. 

We note, however, that we have assumed an unquantified systematic uncertainty in our choice 
of a single prescription for solar convection. While the near surface convection appears not to be 
of great importance to our analysis, we have included no test here of the importance of variant 
convective theories to our results. 

As is the case with stellar modeling in general, the relative importance of element separation 
and mixing is one of our primary uncertainties. Throughout our analysis, we have ignored the effects 
of rotation, except in the decreased efficiency of diffusion encapsulated in the diffusion coefficients. 
In terms of determinations of Tcz, we expect that the most important contributions to be in the form 
of adjustments to the relative efficiencies of mixing and diffusion, which we have shown have an 
impact of Tc^;. For rapidl y rotating ob jects, the rotational splitting of the modes and introduction of 



new modes of oscillation (|Reesell201Cll ) may make mode identification and interpretation challenging. 



Magnetic fields may produces changes in the sound speed near the surface of the star, but we 
expect only a small correction to the sound speed in the deep interior for all but extremely strong 
internal magnetic fields. Again, since the modes of importance have turning points well below the 



photosphere, we expect corrections from magnetic fields to be small. Furthermore, IChaplin et al 



(|2011bl ) finds that highly active stars are less likely to display detectable solar-like oscillations, 
which suggests that the primary role of magnetic fields may be in dictating whether we can detect 
p-modes at all, rather than affecting the acoustic glitch signature itself. 



We have also neglected radiative levitation (see iPinsonneaultl Il997l . for discussion), which 
can selectively levitate some elements relative to others. While our analysis captures the impact 
of global metal diffusion, it does not account for selective levitation of individual elements. In 
particular, this can affect elements such as iron, which contributes substantially to the opacity. 
The effects of radiative levitation are most pronounced in hotter stars with thin surface convection 
zones. The accuracy of the most massive of the models we consider may therefore suffer from our 
neglect of radiative levitation. 

In general, we advocate pairwise comparisons of measured values of Tcz,n for stars which one 
suspects differ significantly in only one way, i.e., testing the mass- relation using two stars of 
very different mass but similar composition and age, or two stars with similar ages and masses 
but different compositions. Obtaining an accurate zeropoint calibration of this relation is currently 
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challenging. For example, variations in Tcz^n on the order of 0.005 can be induced due to numerical 
differences between models with different envelope fitting points. Even in solar models, similar 
numerical uncertain ties due to interpolation can affect the inferred base of the convection zone 
(jBahcall et al.ll2004l ). Physical effects, such as the presence of envelope undershooting, could also 
appear as a zeropoint offset in the relation. Furthermore, as shown in (jBahcall et al.ll2004l ). the 
interpolation of quantities such as radiative opacity tables is uncertain on the 1-3% level near the 
base of the CZ in the Sun. The best approach is therefore to compare pairs of interesting stars, in 
which case zeropoint calibrations will be of less importance. 



Discussion and Conclusion 



We have discussed the factors that affect the location of the base of the convection zone 
thoroughly, but have neglected the second (and probably more commonly discussed) source of 
acoustic glitches in asteroseismic spectra: the helium ionization zone. In the Sun, measurements 
of th e helium ionization zone glitch have constrained the surface helium abundance (IBasu fc Antia 
20041 ). and hopes are high that this will also prove possible in the stellar case. We expect that even 
if the He ionization zone is sensitive to metallicity, the dynamic range of the effect will be much 
too small to make precise metal abundance measurements. Furthermore, the ionization zone lies in 
the outermost layers of the star, and is subject to uncertain surface term corrections, much more 
so than the deeper base of the convection zone. For these reasons a similar analysis on the effects 
of mass and metal content on the helium ionization zone is beyond the scope of this paper, but we 
emphasize that the sensitivity is unknown, and it could yet prove to be an interesting diagnostic. 

Our entire analysis has focused only on main sequence stars, which are inherently fainter, 
and their mode amplitudes are sma ller than the subgaints, which have recently proven to be a ric h 
source of asteroseismic information (iBrandao et al.ll201ll : iMetcalfe et al.ll20ld : iBedding et al.ll2006l ) , 
and are a mong the most comni on stars with detectable solar-like oscillations in the current Kepler 
ensemble (jChaplin et al.ll2011al ) . A preliminary analysis of models evolved onto the subgiant branch 
shows similar trends with composition, with composition effects of similar magnitude. At fixed age, 
Tcz,n now increases with decreasing p and probes models increasingly closer to beginning the ascent 
up the giant branch with deepening convective envelopes. A full analysis of the sensitivity of Tcz,n 
to composition in subgiants is underway. If the sensitivity and theoretical errorbars are similar 
to those on the main sequence, we stand to benefit substantially from extending the analysis to 
subgiants, which have larger mode amplitudes and higher luminosities, which can help to reduce 
observational errors. 



This unique means of measuring the composition promises a host of interesting applications. 
We could, for ex ample, test the tendenc y of planets to be found around hosts of spectroscopically 
high metallicity (jFischer fc Valentill2005l ). Differences between the interiors and atmospheres based 
compositions could help to constrain whether planets are more likely to be found around intrinsi- 
cally metal rich stars, or whether planets themselves tend to enrich the outer layers of their hosts 
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with heavy elements. A simple comparison of spectroscopically and asteroseismically determined 
compositions would in itself be an interesting consistency check, and potentially offer insights into 
the reliability of both methods and physical processes such as element diffusion. The striking sen- 
sitivity of the location of the convection zone to composition even at very low metallicities also 
provides an interesting and relatively rare insight into the interiors of metal-poor stars. These are 
only a handful of the numerous ways in which we can begin to use asteroseismic measurements 
such as this as novel diagnostics of stellar interiors and stellar populations. 

To conclude, we have created a grid of stellar models of different compositions and examined 
the sensitivity of the acoustic depth to the convection zone as a function of composition, mass, and 
our assumptions about the input stellar physics. We make three primary predictions based on the 
analysis of our models: 

1. We predict strong trends in the depth of the convection zone as a function of mass and 
composition. The asteroseismic CZ depth indicator Tcz.n can be different by as much as 
factor of ~ 2 between stars of masses 0.4 and 1.4 Mq. Composition produces changes in 
Tcz,n of order 1% per 0.1 dex in [Z/X]. Tcz,n remains sensitive to the composition even at low 
(~ —1.0) values of [Z/X]. These strong scalings provide a simple test of interiors theory, and 
an absolute abundance measure independent of atmospheric modelling. Furthermore, the 
problem is well posed in Tcz,n — P space, both of which are purely asterosesimic observables. 

2. Reasonable estimates of theoretical and observational uncertainties suggest that not only is 
'Tcz,n sensitive to the composition, but that the uncertainties in the relationship are small. On 
average, we expect to be able to measure absolute abundances to 0.15- 0.3 dex for solar- like 
stars at 5 Gyr given the assumed observational and theoretical uncertainties. 

3. Finally, the measurement of the depth to the convection zone has potential diagnostic power 
as a means of probing theoretical uncertainties. In particular, we have addressed the manner 
in which one would use Tcz.n to test for rotational mixing in Li dip stars, and to test for 
differences in the relative element abundances in an ensemble of targets. 

Measurements of Tcz,n have the potential to both constrain interiors theory in terms of the 
balance between diffusion and mixing, element abundance patterns, and the basic prediction of a 
strongly mass dependent CZ depth. The technique also and offers a unique, absolute abundance 
measure, which is inherently useful in the study of the chemical enrichment of the galaxy, and 
benchmark for comparison to stellar atmosphere derived abundances. This is a powerful tool that 
can help us to precisely measure stellar parameters and test the physics of stellar interiors. 
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Table 1. Theoretical Errorbar Model Grids 



Grid name 


a 


X 


Y 


Z 


Description 


standard 


1.93271 


0.710040 


0.018338 


0.271622 


standard physics, Section 12.21 


mixing 


1.80975 


0.719426 


0.016468 


0.264107 


no diffusion 


eos (scv) 


1.87640 


0.709470 


0.018455 


0.272075 


SCV/Yale EOS 


eos (yale) 


1.90600 


0.716180 


0.018665 


0.265155 


Yale EOS only 


overshoot 


1.91560 


0.707060 


0.018299 


0.274641 


Convective core overshoot 


mixture 


1.90600 


0.720880 


0.014760 


0.264360 


AsDlund et al. f2009') solar mixture 


nuclear 


1.92500 


0.709110 


0.018257 


0.272633 


nuclear cross sections Sll, S33, S34, S114 +4cr 


opacity 


1.92250 


0.709415 


0.018338 


0.272247 


OPAL/ Alexander opacities 


alpha 


many 


0.710040 


0.018338 


0.271622 


a values, chosen to produce AR = +0.01 


boundary 


1.82100 


0.710030 


0.018334 


0.271636 


Grey atmosphere boundary condition 


diffusion 


1.95000 


0.708850 


0.018597 


0.272553 


diffusion coefficients altered by 15% 



M (MJ 
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M (Ma 
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Fig. 1. — Left panel: The physical depth of the convection zone, normalized by the radius of the 
star, as a function of mass. Middle panel: The acoustic depth to the convection zone, normalized 
by the acoustic depth from surface to center of the star. Right panel: The normalized acoustic 
depth as a function of effective temperature. The dotted curve is for solar composition models at 
1.0 Gyr, the solid for 5.0 Gyr, and dashed for 10.0 Gyr. All models are on the main sequence, with 
the central hydrogen fraction Xc > 0.0002. The top axis in the center and right panels gives the 
mass for objects at 5.0 Gyr (solid curve). The corresponding discussion can be found in Section 

EH 



Table 2. Model Grid 



Mass Yi 


[Z/X], 


[Z/X] 


Age 


log(^) log(^) T,ff log(^) R,,/R Tc. 




Tcz,n 




(initial) 


(surf) 


(Gyr) 


(K) (s) 


(s) 





0.400 


0.240 


-1.2 


-1.20e-f00 


0.500 


-1.611 


-0.4313 


3754 


0.8961 


0.6090 


839 


1237 


0.6791 


0.400 


0.240 


-1.0 


-l.OOe+00 


0.500 


-1.630 


-0.4306 


3711 


0.8940 


0.6032 


847 


1240 


0.6834 


0.400 


0.240 


-0.80 


-8.01e-01 


0.500 


-1.652 


-0.4308 


3664 


0.8946 


0.5926 


854 


1238 


0.6908 


0.400 


0.240 


-0.60 


-6.01e-01 


0.500 


-1.676 


-0.4328 


3623 


0.9005 


0.5736 


865 


1229 


0.7046 


0.400 


0.240 


-0.40 


-4.01e-01 


0.500 


-1.697 


-0.4361 


3592 


0.9105 


0.5444 


880 


1214 


0.7254 



''Table [2] is available in its entirety at www.astronomy.ohio-state.edu/CZdepth/model_grid.txt. A portion is shown here for 
guidance regarding its form and content. The columns are as follows: 1) mass in solar masses, 2) initial helium abundance, 
3) initial model [Z/X] referenced to the initial Z/X of a calibrated solar model, 3) the surface [Z/X] at a given age, 4) age, 
5) log luminosity in solar luminosities, 6) log radius in solar radii, 7) effective temperature, 8) log mean density given by 
log{MR~^ /MqRq^), 9) fractional radius of the CZ, 10) acoustic depth to the CZ, 11) acoustic crossing time, 12) normalized 
acoustic depth. 
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Te« (K) 

Fig. 2. — The variation of the normaUzed acoustic depth of the base of the convection zone as a 
function of composition and T^ff. Sohd hnes represent models of fixed initial [Z/X] (reference to 
the initial solar abundance) from red/top most (+0.6) to purple/bottom most ( — 1.2) spaced every 
0.2 dex. The composition dependence of the acoustic depth is plotted for three representative ages, 
and dotted lines are plotted in gray for constant mass, in solar units. All models are stars on the 
main sequence, with Xcore ^ 0.0002 and have an initial, solar-calibrated helium of 1^ = 0.271. 
Representative observational error bars on the c[ua<ntities Tcz,n 

and Tef / are shown in the left most 

panel. See Section [3r2l for discussion. 
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Fig. 3. — The variation of the absolute acoustic depth (in seconds) of base of the convection zone as 
a function of composition and T(,ff. Sohd hues represent models of fixed initial [Z/X] from red/top 
most (+0.6) to purple/bottom most (—1.2) spaced every 0.2 dex. The composition dependence of 
the acoustic depth is plotted for three representative ages, and dotted lines are plotted in gray for 
constant mass, in solar units. All models are stars on the main sequence, with Xcore > 0.0002 and 
have an initial, solar-calibrated helium of 1^ = 0.271. Representative observational error bars on 
the quantities Tcz and Teff are shown in the left most panel. See Section [32] for discussion. 
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Fig. 4. — The variation of the normahzed acoustic depth of the base of the convection zone as 
a function of composition and p. Sohd hnes represent models of fixed initial [Z/X] from red/top 
most (+0.6) to purple/bottom most (—1.2) spaced every 0.2 dex. The composition dependence of 
the acoustic depth is plotted for three representative ages, and dotted lines are plotted in gray for 
constant mass, in solar units. All models are stars on the main sequence, with Xcore > 0.0002 and 
have an initial, solar-calibrated helium of 1^ = 0.271. Representative observational error bars on 
the quantities Tcz,n and p are shown in the left most panel. See Section 13.21 for discussion. 
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Fig. 5. — The variation of the absolute acoustic depth (in seconds) of base of the convection zone 
as a function of composition and p. Sohd hues represent models of fixed initial [Z /X] from red /top 
most (+0.6) to purple/bottom most (—1.2) spaced every 0.2 dex. The composition dependence of 
the acoustic depth is plotted for three representative ages, and dotted lines are plotted in gray for 
constant mass, in solar units. All models are stars on the main sequence, with Xcore > 0.0002 and 
have an initial, solar-calibrated helium of 1^ = 0.271. Representative observational error bars on 
the quantities Tcz and p are shown in the left most panel. See Section 13.21 for discussion. 
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Fig. 6. — Left panel: The fractional difference in Tcz,n between the surface abundance and initial 
composition at constant [Z /X] . This difference arises because of gravitational settling of heavy ele- 
ments, and would manifest itself as a T^ff dependence of the surface [Z/H] in in mono-composition 
sample. Center panel: Lines of iso-composition in surface abundance, compared to a reference 
model at solar surface abundance. Left panel: surface abundance as a function of p for models of 
solar composition. The solid line denotes models at 1 Gyr, the dotted at 5 Gyr, and the dashed at 
10 Gyr. See Section 13.21 for discussion. 




Fig. 7. — The fractional difference in Tcz,n between models with an initial solar helium abundance 
{Y = 0.271) and Y = 0.24 (sohd), Y = 0.26 (dotted) and Y = 0.28 (dashed) at constant initial 
Z/X (solar) for models of 1.0 Gyr. ^Tcz,n/'Tcz,n > represents models in which the convection zone 
is deeper than in the standard, solar Y case. See Section 13.21 for discussion. 
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Fig. 8. — The fractional uncertainty in the measurement of Tcz^n for a 5 Gyr old, solar composition 
star from our each class of uncertainties, from left to right: random, systematic, and observa- 
tional. See Section \T3\ for a discussion of the uncertainty calculations. The sources of uncertainty 
are coded as follows: left panel: solid line: total, dotted: diffusion coefficients, dashed: nuclear 
reaction rates. Center panel: solid- total, dotted- EOS, dashed-overshoot, dot dashed-opacity, 
double-dot dashed-boundary conditions. Right panel: solid-total, dotted-Y, short dashed- age, 
dot dashed-measurement of p, triple-dot dashed- measurement of Tcz,n, long dashed-mass-radius 
relation, through changes in a. 
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Fig. 9. — The uncertainty in the measurement of the metahicity of a solar composition, 5 Gyr 
old star using Tcz,n as a composition indicator. The dashed (blue) line represents the contribution 
from asteroseismic uncertainties on the observables, the dashed (purple) from observational uncer- 
tainties on the mass-radius relation, Y, and age, the triple-dot-dashed (orange) from systematic 
uncertainties, long-dashed (red) the zeropoint uncertainties (mixture + mixing), dot-dashed (green) 
the random uncertainties, and solid (black) the combined systematic, random and observational 
uncertainties. See Section 13.31 for discussion. 
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Fig. 10. — The fractional difference in Tcz^n (the normalized acoustic depth to the base of the 
convection zone) between models of different compositions as a function of T(,ff. Each solid line 
represents the fractional difference in Tcz,n between a given [Z/X] (—1.2, —0.8, —0.2, 0.2, 0.6, marked 
for reference in the right panel) and Tcz,n for a solar composition model. The gray shaded region 
represents observational and theoretical errors on Tcz,n, both described in detail in Sections 12.31 
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Fig. 11. — The fractional difference in Tcz,n (the normalized acoustic depth to the base of the 
convection zone) between models of different compositions as a function of p. Each solid line 
represents the fractional difference in Tcz,n between a given [Z/X] (—1.2, —0.8, —0.2, 0.2, 0.6, 
marked for reference in the right panel) and Tcz^n for a solar composition model. The gray shaded 
region represents observational and theoretical errors on Tcz,n, both described in detail in Sections 
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Fig. 12. — Fractional difference in Tcz,n between a set of models with standard physics, and a set 
where diffusion has been eliminated entirely to mimic efficient mixing. The solid (blue) line is for 
models at 1 Gyr, dashed (green) at 5 Gyr, and dot-dashed (red) at 10.0 Gyr. Note that over the 
temperature range in which the Li dip is observed, the difference is of order 6%, such that over a 
very narrow temperature range we expect to see and aburpt change in the location of the the base 
of the CZ. For the 5.0 and 10.0 Gyr curves, all stars in the Li dip temperature range have already 
evolved off of the MS. See Section [3.4.11 for further discussion 



